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In a paper read before the Society on June 11, Sir William Thomson expressed a 
doubt as to the general truth of the Maxwell-Boltzmann doctrine concerning the 
distribution of energy among a great number of mutually acting bodies, and suggested 
that certain test cases should be investigated. The test that he proposed on that 
occasion was a number of hollow elastic spheres, each of mass M, and each containing 
a smaller elastic sphere of mass m, free to move within a larger one. This pair he 
calls a doublet. This case is within the general proof of the doctrine given below. 
It is, however, I think amenable to a simpler treatment, which has been applied to 
the case of elastic spheres external to one another. 

1. Every doublet has a centre of inertia of the sphere M and its imprisoned m. 
Let Y be the velocity of that centre of inertia, It the relative velocity of M and m. 
If V and It be given in magnitude, It given in direction, Y may have any direction, 
and in Maxwell’s distribution, for given direction of It, all directions of V are 
equally probable. Conversely, if, whatever be the values of Y and It, for given 
direction of It all directions of Y are equally probable, Maxwell’s law prevails. 
Now consider a very great number of doublets, all having their relative velocity and 
the velocity of centre of inertia within limits R, It -j- dR, and Y,. V -f dV. Consider 
them before and after collisions between M and m. Nothing is changed by collision 
except the direction of It, and that change of direction is independent of the direction 
of Y. Therefore after collision for given direction of It all directions of V are equally 
probable, and therefore Maxwell’s distribution prevails after as well as before col¬ 
lision, and is therefore not affected by collisions. 

To proceed to more general cases. 

2. The characteristic property of collisions of conventional elastic bodies is that 
with continuous variation of the coordinates, and without variation of the kinetic 
energy, the velocities at a certain instant change discontinuously. The general 
treatment adapted to systems of this kind is as follows:—Let there be a system 

defined by n coordinates p l . . . p n , the corresponding velocities being p 1 . . . p„, and 
the generalised components of momentum q± . . . q„. At a certain instant a collision, 

fie., discontinuity of p 1 . . . p' n occurs. After the collision, let the velocities and 
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momenta be denoted by p\ . . . p' n , q\ . . . q n . Let U be the potential. In the con¬ 
figuration p x . . . p n there are n — 1 independent linear functions of the n forces, 
dTJ dJJ 


dPi' ' ' 
such that 


dp, 


, each equal to zero. We might then find n new coordinates < 


. • Cm 


(HJ 

dCy 


0 . . 


dc n . 


0, and therefore if Sj, S 2 . . . S„ be the components 

dS,, 


p . dSi dT 0 ^ dSi 

or momentum corresponding to c x . . . c n , -37 = — &c., and —- 


-1 


dt d.c Y ’ dt ’ ' ’ dt 


are finite or 


zero. In the limit when the remaining force becomes infinite and acts for an infinitely 

short time t, f ^ dt = 0 or S'j — Sj = 0 ... S',,-! — S,,-! = 0, a,nd restoring the 
J o dt x 

original coordinates. 


«i iPi ~ P\) + \ (Pz~p\) + ... + h (_ p n -p „) = 0‘ 
a z(Pi -P\) + • • • = 0 




= 0 


(A), 


«*- (Pi - Pi) + • • • 
in which the coefficients are functions of the coordinates. 

And therefore n — 1 li near functions of the p’s are unaltered by the collision, namely 


Therefore 


%p'q = tpq. 


2(?+ 9’) (P ~ P) = 0 


a iPi + hP» + • • • + h Pn = a iP\ + hPz + • • • + hPn = suppose"! 

a zPi + l } dh + ■ • • + h Pn = c hPi + P\ + • • • + hp'n — 

a n-lPl 4" b n -\Pn 4" • • • 4" h-lPn — C^n-lPl 4" ^n~lP'z 4~• • • 4" & n-\P n— S*-! 

Again, since the kinetic energy is unchanged, 

tpq= Zp'q', 

and by the properties of generalised coordinates 


(B). 


■ (C). 


The last equation forms, with the n — 1 equations (A), a system of n equations, 

all linear as regards p x — p\, &c. Since p x — p\, &c., are not all zero, we must 
equate the determinant of the system (A) and (C) to zero. That gives us a linear 
equation in q x + q\, q % + q'%, &c. We can now substitute for the q’s their values in 

terms of p x . . . p n , and so obtain a linear equation 









Downloaded from rsta.royalsocietypublishing.org on July 5, 2012 


MR. S. H. BURBURY ON THE COLLISION OE ELASTIC BODIES. 409 

a (i>i + Pi) + P(Pa+P'a)> &c - = 0, 
or 

a Pi + &p» + ■ • • = R = — ( a p\ + ftp* • • •)» 

where a, /3, &c. are functions of the coordinates and constants. 

That is, we have in all n — 1 linear functions of the velocities, namely, S x ... S«_u 
which remain unaltered by the collision, and one other linear function, It, which 
remains unaltered in value, but changes sign. That must be the case on every 
collision of elastic bodies. 

3. The kinetic energy may be expressed in terms of the n variables S x . . . S w _ a , 

and It, in lieu of the n velocities yq. .. p n , and since it is not altered by the collision, 
which changes the sign of It, leaving S x ... S„_j unaltered, it must be of the form 

2E =/(Sj... S„_i) + XR 2 , 

in which/(Bj.. . S»_j) is a quadratic function of S x .. . S»_i with coefficients functions 
of the coordinates and constants of the system, and X is a function of the coordinates 
and constants. 

4. The system, after collision, has velocities p\, &c., which we will call the 
second state. We may conceive a system with the same coordinates having 

velocities — p\, — p\, &c., and call this the second state with reversed velocities. In 
this state, Sj.. . S«_! will have opposite signs to those they have in the first state, 
and R has the same sign as in the first state. The system retraces its course, and a 

collision occurs changing — p\ into — yq, &c., leaving Si. . . unaltered, and 
changing R into — R. 

5. To define a collision, we may suppose that a certain function, xp, of the coordinates 
and constants is prevented by the physical conditions of the system from becoming 
positive. When xp becomes zero, dxp/dt from being positive becomes discontinuously 
negative, and a collision is said to take place. We may take xp for one of our 

generalised coordinates in lieu of p„, and ip, or dxp/dt, for the corresponding component of 
velocity. The kinetic energy is a function of p x - • • Pn-i, 4>, and we may express it as 
a function of S 2 .. . S„_,, and ip, where S x . . . S B _! are the constants found above. 
Since the kinetic energy is not altered by the discontinuous change in ip, whatever the 

values of S t . . . S«_ l5 it must be of the form f(S 1 ... + i|Xt/d. That is ip is 

reversed in sign, but unaltered in magnitude by the collision, and is, therefore, equal 
or proportional to R found above. 

6. What we have proved for a system is of course true if for system we write pair 
of systems. For instance, let there be two sets of systems: (1) systems M defined by 

coordinates yq . . . p r , and velocities yq . . . p r , and (2) systems to defined by coordi 

nates and velocities yq +1 . . . p n and y> r + 1 . . , p n . If xp is a function of p y ... p n such 
MDOCCXOII.-—A. 3 G 
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that when xj/ = 0 a collision— i.e., a discontinuity of the velocities—occurs, we may 
treat the pair of systems M, m in all respects as a single system within the preceding 
investigation. 

7 . All those systems, or pairs of systems, for which at any instant \p lies between 
zero and — ( d\\>jdt ) St, cl\\>jdt being positive, will undergo collision within the time St 
after that instant. We may, therefore, take dxjf/dt or It as measuring the frequency 

of collisions for given values of p u &c. 

8. From the linear equations (B) above given we can find any of the velocities, for 
instance, jq, as a linear function of S x . . . S»_ x R, and p\ will be the same function 
with — R written for R. Hence jq 2 — p\ 2 = 4R2/i,S where 

YjuS = /qSj + /a 2 S 2 + 

and the pa are functions of the coordinates and constants. 

Also (yq 3 — p\ 2 ) R = 4R' 2 2juS. 

Row, without altering E or R, or the coordinates, let us make S x . . . S„_i pass 
through the whole range of values consistent with 

2E = \ 3 R 2 +/(S x . . . S„_0.(E). 

Als o let (j> (S x . . . S„_j) dS L . . . cZS„_x be the number of systems for which these 
variables lie between the limits 

S x and Sj + dS, 

S,_ x and S M _i ~h dS„_x, 

E and R and the coordinates being constant. 

Then 

ff... (k - p'i ) R • $ (s x ... s H _x) ds x ... ds M _x 

= 4R 3 jj. . . <£ (S x . . . S*_x) tp S dS x . . . dS„_„ 

the integrations being over all values consistent with (E). 

Now, in the Maxwell-Boltzmann distribution <f> (S x . . . S»_x) is a function of 
the kinetic energy only, and is therefore constant throughout this integration. There¬ 
fore the right-hand member of the last equation is zero, because for any given set of 
values of Sj . . . S„_i satisfying (E), there is a corresponding set with the opposite 
signs also satisfying it. Therefore 

j|. .. R<£ (Sj. . . S*_x) (p* - p'p) dSj . . . d8n-i = 0, 

and, therefore, the average value of jq 2 — p x 2 for all collisions given E and R is zero. 
The same is true for yq 3 — p &c. 
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And, therefore, since E and R are arbitrary, p 2 , &c., are not altered by collision 
at all, that is, the Maxwell-Boltzmann distribution, given existing, is not altered 
by collisions. 

The above proof also shows that p 3 — p' 2 is zero without the factor 11, that is, the 
average value for all systems is zero, as well as for all collisions; and in proving that 

p 3 = p' 2 , it does not matter whether we introduce the factor R or not. 

9. We will give certain examples of the functions S x . . . S„_ x R. 

I. Elastic spheres of masses m and M respectively. Here a colliding pair, which 
corresponds to a system in the general treatment, has six degrees of freedom. There 
should, therefore, be five linear functions of the velocity unaltered by collision. They 
are x, y, X, Y, the tangential components of velocity at the instant of collision, and 

mu + MU = (M + m) Y = mu' + MU', 

where u, U are the normal components. 

We have from the last equation 

m(u~u') + M(U-U') = 0, 

and by the equation of energy 

m{u 3 - u' 2 ) + M (U 3 - U /3 ) = 0, 

whence 

u + u' = U +U', 
or 

u-TJ-- (u - U') = R. 

II. The system consists of a sphere of mass m colliding with a spheroid of mass M. 
It is assumed that the spheroid will acquire no rotation about its axis of figure, but 
may have rotation about any other principal axis. It has then five degrees of 
freedom, and the system of sphere and spheroid has eight. 

We require, then, seven linear functions of the velocities to be invariable. 

z 



Let 0 be the centre of the spheroid, OZ in the plane of the paper, its axis of 
figure, P the point of collision, PN normal at P, ON perpendicular to PN, and 
ON — c, A the moment of inertia of the spheroid round an axis through O perpen¬ 
dicular to OZ. Then our seven constants are 

3 G 2 
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(1), (2), (3), (4). Four tangential components of velocity. 

(5). The angular velocity 6 of the spheroid round an axis perpendicular to OZ in 
the plane of the paper. 

And if u, U be normal components of velocity, and co the angular velocity round an 
axis through 0 perpendicular to the plane of the paper, the following two, viz. :— 

By conservation of momentum, 

(6.) mu + MU = (M + to) Y = mu' + MU'. 

By conservation of moment of momentum round the axis through 0 perpendicular 
to the plane of the paper, 

(7.) - McU + Aco = S = - McU' + Aw'. 

We then form the equations 

m(w-0 + M(U-U') = 0, 

- Mc(U - U') + A(w~ co') = 0, 

and by conservation of energy, 

m (u + v!) (u - n) + M(U + U')(U — U') + A (a> + o>') (<» - «') = 0 ; ■ 

and equating the determinant to zero we obtain, neglecting common factors, 

u — U — cco — p, 

/ TTY / 

u — U —— Co) — p P 

and, using V, S, and p for . . . S„_i and II in the general equations, we see that if 
the Maxwell-Boltzmann distribution of energy exist, it is not disturbed by colli¬ 
sions between spheres and spheroids. 

III. Professor Burnside’s problem (see his paper, Roy. Soc. Edinburgh, July 18, 
1887). He supposes a number of similar and equal spheres, each of unit mass, but 
each sphere, instead of being homogeneous, has its centre of inertia at a distance c 
from its centre, c being supposed very small compared with the radius. The principal 
moments of inertia are for each sphere A, B, C, and the direction cosines of c referred 
to the principal axes through the centre of inertia are for each sphere the same, viz., 
a , /3, y. The direction cosines of the line of centres at impact referred to the principal 
axes are for one sphere L, M, N, and for the other l, m, n. 

Further the normal velocities are U, u, and the angular velocities round the 
principal axes are XR, fl a , fl 3 , co v co z , w 3 for the two spheres respectively. 

Finally we write, 

NyS — My — P n/3 — my = p. 

Ly — Na — Q ly — no. = q. 

Ma — L/3 = R ma. — = r 
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The system of two spheres has twelve degrees of freedom. We require, there¬ 
fore, eleven linear functions of the velocities to be invariable. They are as follows, 
viz., four components of velocity in the common tangent plane, x, y, X, y, and seven 
others, viz., 

u + U = V = u' + U', 
epu + Aft)! = s x = epu + Aft)'!, 
cqu -p Bft) 3 = s. 2 = cqu -p Bco\, 
cru + 0 ft) 3 — s 3 — cru + C co' s , 
cPU - Aip = Si = CPU' - Aa'j, 
cQU - Bn 3 = S 3 = cQU' - Bft' 3 , 
cRU - Cfl 3 = S 3 = cRU' - Cft' 3 . 


As before, we form the equations, 

u — u' + U — U' = 0, 

cp ( u — u) + A (ft)j — w i) = 0, 

cq (u — u) + B(c«) 2 — ft/ 2 ) = 0, 

cr (u — u') + 0 (ft) 3 — ft.'s) = 0, 

cP(U-U') — A (il x — il'J = 0, 

cQ(U-U') — B (0 3 — il' 3 ) =0, 


cB (U - O') - C (il 3 - il' 3 ) = 0, 

and by the conservation of energy, 

(u + u) (u - u') + (U + U') (U - U') + A (ft)! + ft)'i) (ft)i - ft/j) + &c. = 0 
and equating the determinant to zero, 

u — U — c (jpoi i -p q(o% -j- ^ft) 3 -p Pili ”1“ Q ^2 "p B;il 3 ) — p, 
it — U' — c (pco'i + qco'z + r&)' 3 -p Pil'i + Qil'a + Hil' 3 ) = — p 

We can now substitute V, .q, q, q, S x , S 2 , S 3 , p for S x . . . S„_iR in the general 
equation, and we obtain, as before, the result that the Maxwell-Boltzmann distri¬ 
bution, given existing, is not affected by collisions. 

10. Professor Burnside obtains the same eight equations as above given, and I 
acknowledge my obligation to him, but he originally deduced the result that the 
energy of rotation is twice the energy of translation, instead of equal to it, as, 
according to the theory, it should be. He has since seen reason to change his views 
with regard to this problem. 
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The equations given by Professor Burnside can easily be modified so as to meet the 
case of elastic bodies of any shape. 

11. It would not be difficult to extend the method of (8) and show that the 
Maxwell-Boltzmann distribution is a necessary, as well as a sufficient condition for 
stationary motion. But that is more completely done by following or extending 
Boltzmann’s method. 

Let there be a set of systems which we will call systems M, whose co-ordinates and 
velocities are p x ■ ■ ■ p r and p x . . . p,.. 

Let the number of such systems for which at any instant the co-ordinates lie 
between the limits 


p x and p x + dp x 

>r + dp, J 


p r and p r 


(A), 


and the velocities between the limits 


be 

or, shortly, 


p x and p x + dp x 
'p r and p r + dp,. 

F (p x . . . p r p x . . . p } ) dp x . . . dp r dp x . . . dp r , 

¥dp y . . . dp r dp x . . . dp r . 


■ (A'), 


Let there be another set of systems, which we will call systems to, whose co¬ 
ordinates are p r+1 . . . p n and velocities p r+x . . . p n . 

And let the number of systems to for which at any instant the coordinates lie 
between the limits 


Pr + i and p r+1 + dp r + l " 

.>- 

p n and p n + dp H 


(B), 


and the velocities between the limits 


be 


Pr+i and p> r + j -f dp r + l 

p„ and p n + dp„ 

f(Pr +1 • • • PnPr + l • • ■ pn) dp r + 1 . . . dp„ dp r+ , . . . dp„ 


• (B'), 
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or, shortly, 


/• dp. 


+1 


dp n dp r 


+ i 


dp n 


415 


The number of pairs of systems, each consisting of one system from each set, whose 
coordinates and velocities at any instant lie within the above limits, is 

dp y . . . dp r dp r+l . . . dp n F f 

Now let i(i be a function of the coordinates which cannot become positive, and 
such that when xp = 0 the velocities change discontinuously, and a collision occurs. 
We may use t p for one of our coordinates, expressing p n in terms of p } . . . p „_, 

and xp, and in like manner p n in terms of P\ • • • p n -\ and xp. All those pairs for which 

at the given instant xp lies between zero and — xpSt, xp being positive, will undergo 
collision within the time St after that instant. And so the number of such collisions 

in unit time is, writing It for xp, 

dpy. . . dp n _ 1 dpy . . . dp n _j c£R.. F/R. 

Here F is a function of p x . . . p r only, but/*, by virtue of the elimination of p n and 
p n , is now a function p x . . . p n -i and p x . . . p n _ x R. 

Each of these collisions changes p x into P\, &c., that is, changes the system from 
the first into the second state. The number in unit time of collisions, which with the 
same coordinates change the system from the second, state with reversed velocities 
into the first state with reversed velocities, is 

dpy . . . dp n _y dp\ . . . dp' H _ 1 cHtF'/Tt, 

■A 

in which F', f are the same functions of p\, &c., as F,/are of py, &c. 

By a known theorem 

dp'y . . . dp'xi _! = dpy . . . dp„_y, 

and in Maxwell’s distribution F f = F/'. And so the number of reverse collisions 
is in that distribution equal to the number of direct collisions. And this insures the 
permanence of the distribution. It is assumed that there are always as many systems 
with any given set of velocities as with those velocities reversed. And so we speak 
of the second state with reversed velocities as equivalent to the second state. 

12. If F/' =f= 1/ the number of reverse collisions is not equal to the number of 
direct collisions. And, therefore, more (or fewer) pairs of systems pass out of the 
first state into the second than vice versa. In that case the number of systems M, 
whose coordinates and velocities are within the limits A, A', is increased by collisions 
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with systems in, whose coordinates and velocities lie within the limits B, B', by the 
quantity 

dpi • • • dpn -1 dpi ■ ■ • dp u _ l c/R (F f — F '/) R per unit of time, 

and is increased per unit of time by collisions with systems to for all values of 
p,. + 1 . . . p„ by the quantity 

dp 1 . . . dp r dp x . . . dp r |j" . . . (Ff — Ff) R dp r + 1 . . . dp n _ x d.Jx, 


all values of p r+ i, &c„, being included in the integration. 

We will now assume (see p. 418, post) that the velocities of M and to are not on 
average altered except by collision between M and to. And so the above-mentioned 
increments are the only increments by which the class of M systems within the limits 
A, A' is affected. 

In that case 


ilP 

dt 


1 • . r r 

dp l . . . dp r = dp l . . . dp,. . . . (F/' — F/) R dp r+l . . . dp tl _ x dR, 
and, therefore, 

||... log F dp x . . . dp r (over all values of p x . . . p,) 

= ||. . . (F'f - Ff) R log Fdp 1 . . . dp n _ x dK 


By symmetry, as the right-hand member includes all possible collision between 
M and to, 


ff • ■ • + • • dp, 

= ff •, • (F'f' - v/) E log/#,... dR, 

and, therefore, 

If • ■ • if I °g I ’ ti Pi • • • dp, + jf . . . " log/#, + 1 .. . dp. 
= ff . . . (F'f - Ff) E log (Ff) dp,... dp,., dR. 


By symmetry, as we may interchange the accents, 

ff • - • fr log F dpi • • • dpr + ff ’ ' ' % l ° g f dp r + \ ■ . . dp„ 

= HI ...(Ff- Ff) R log (Ff) dp y . . . dpfi dll, 
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and, therefore, 

= i fff ... (vy ' - I f) E log |£ dp,... dp,., m. 

Now let 


H = fff • • • F 0°g F - 0#, • • • dPr + fff • • ../(log/ — 1 )dp r+1 . . . dp„ 


and, therefore, 


dl I 

dt 


= fff • • • ]§ lo g F d Pi ■ ■ ■ dpr + fff • • • f log/ dp , + 1 . . . dp„ 
= | fff. .. (F /' - F/) R log dp,. . . dp,., dU, 


as we have seen. 

Now, this expression is necessarily negative, unless F f — F f whenever the pair of 
systems having coordinates and velocities p l . . . p n can pass by collision, and, there¬ 
fore, with unchanged kinetic energy, into the state in which they are p x . . . p',„ that 
is, unless the Maxwell-Boltzmann distribution exist, and is then zero. H therefore 
tends to a minimum which it reaches when F/ = F 'f. 

We will now make 

H - H x + K 


where IIj is the minimum value assumed by H when F '/' = iy! Then dH./dt = 0, 
and dK/dt = dH/dt, and is always negative. We may define the function K to be the 
disturbance, and (l/K) (dK/dt) to be the rate of subsidence of the disturbance by 
collision. In certain cases, we can calculate the rate. 

13. We have assumed that f varies only as the result of collisions. That is, if 
of ft denote the time variation of f due to causes other than collisions, and fill ft be 
formed from dfft as dH/dt from df /dt, then dH/dt = 0, on average. It is worth 
while to consider on what condition this may be safely assumed. 

Let 


dt 



log/c/p! ... dp„. 


As we are dealing with rigid elastic bodies under the action of no forces, we may 
treat f as a function of three translation velocities, and three angular velocities, 
w 1 , iv 2 , iv z , about three principal axes of the solids. Let A, B, C, be the principal 
moments of inertia. Evidently, there being no forces, the translation velocities cannot 
vary except as the result of collisions. But for each solid, w x , iv z , w % may vary, the 
law of their variations being Euler’s equations. We may, therefore, in calculating 
0H ft treat f as a function of io u iv 2 , iv ?> only. Then 
MDCCCXCII.—A. 3 H 
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df df dw l 

dt dw l dt 


+ 


df dw 2 
dwc, dt 


+ 


df dvi, 
dw s dt 


B - C df 
A du\ 


w 2 tv 3 


+ 


C - A df 
B dw % 


, A - B df 

+.o.. w i 


with two other corresponding terms, the limits being in each dr 00 • Now, with these 
limits. 


\\\faf lo g/- w * v h ■ dw i dv h dl h 

= J] dw 2 dw $ W 2 V f {/log/. 


■/log/ 




Now, we may assume f—0 and f log f = 0, when any one of the three variables is 
infinite, whether positive or negative. And this assumption is sufficient to justify 
the statement 3H/9i = 0. 

14. It is possible to calculate, in a simple case, the rate at which a disturbance subsides 
by collisions. For example, two sets of elastic spheres, N of mass M, and n of mass 
m, in unit of volume. In the normal state, the number in unit of volume, whose 
velocities are represented by lines drawn from an origin to points within an element 
of volume U 2 sin 0 d0 df d\J is for the M spheres 

N J e - ; ' MD2 U 2 sin 0 d0 df dU, 


where U, 0, f are usual spherical coordinates ; or, let us say, 


F (U) = N 



AMU-* 


Similarly, for the m spheres, 

/(w) = n 



e - him i 2 


expresses the law of distribution of velocities in the undisturbed state. We will 
write F and f for these expressions. 

We will now suppose there is a small disturbance consisting in h having different 
values for the two sets. Let h be written h (1 + F>) for the M spheres, and h (l + d) 
for the m spheres. We shall neglect third and higher powers of D, d. Then F 
becomes in the disturbed state 




Downloaded from rsta.royalsocietypublishing.org on July 5, 2012 


MR. S. H. BURBURY ON THE COLLISION OF ELASTIC BODIES. 

that is, 

N (~j (1 + D)s e - ; ' MU2 $ i _ ADMU a + terms in I) 3 }. 

(It will appear that terms in D 3 , &c., are not required.) 

Similarly /becomes 


n 


(v)’ ( L + # e~ hmu2 (l 


— dhmu 3 + d 3 &c.). 


419 


We will further suppose that the disturbance is introduced without changing the 
total energy. That gives the relation 


N, 


-|- n 


2h (1 + D) 1 2h (1 + D) 2 h 




or 


N 


+ 


1 + D 1 1 + d 


= N -f- n. 


15. The disturbance will subside by collisions between M and m. And we will 
treat of the case in which it subsides in such manner that the above values of F and/ 
apply at every instant with the values that D and d at that instant have. Such a 
mode of subsidence is possible, at all events if our equations lead (as they do) to 
a relation of the form (1/K) (dK/dt) = constant. 

Let us then form the function 


h= nr r. f (i og f -1) u 3 S m «dp ^u, 

J 0 j 0 j 0 


where U, a, and /3 are usual spherical coordinates, 


+ 


G JOJO 


/(log/ — l) id sin a da d/3 du, 


where F and/have the values above given. 
That is 


or 


because 


H = f log (1 + D)|[[f.U 3 sin a da d/3 dU 

+ t i°g‘ (1 + d) HI f.u~ sin a da d/3 du 
-j- terms independent of D and d, 

H = N f log (1 + D) + n flog (1 + d), 

+ terms independent of D and d, 

|j[FU 3 sin a da d/3 dU = N 

|| f f u '~ ®i n a do. dp du = n. 

3 H 2 
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The terms independent of D and cl are equivalent to H 1; the minimum value of H 
when there is no disturbance. 

And so 

K = H — Hj = f N log(l + D) + fnlog (1 + d). 


That is 


Now since 


K = fn^-f) + fN(D-^y 


D 2 ' 


and 


N i n AT I 

l + D + n r n ’ 

nd 


D 


K = # n [ cl — 




N. 


N | (N + n) d 


nd 




N + (N + n) d 4 ’ {N + (N + n) d}°~ 




f . (1 - - ~- <l 


n~d~ 


4 {N + (N + n) d}* 




N 


— 3 


N 


(N + n) d 2 . 


In order to find dK/dt we will transform our coordinates thus : Let V denote the 
velocity of the centre of inertia of a pair of spheres M and to, p their relative velocity, 
6 the angle between V and p. Then 


U 3 = V 3 + 
v? = Y 2 + 


m 


P + 


2m 


Yp cos 0 


pl— , , , Yp cos 0 
M + m' ) M + m 


M + m ' J M + m 
M \s 2M 


and 


and so 


MU 3 + mu 2 = (M + to) V 3 + 


Mm 
M + m 


F/= N«(v) # (vj(l + D )*(l + 4^- A { (M + w)V2 + smr/i (1 - DAMU 2 - dhmv?) 


Yf = Nw (— Y ( Am Y (1 + L>) f (1 + dy e- A (< M + -> V3+ imr/} (1 - D/iMU' 2 - dhmu 2 ) 


where U', u' are the values of U, u after collision. 
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Therefore 


F f- Ff 


Nn (~J (v)* ( l + D)*(l + d)» + 




{DAM(U' a - U 2 ) + dhm (id -id)} 


Now, if 9 ' he what 9 becomes after collision, 


U' 3 - IP = 2Y p . 


m 


M + m 


(cos ff — cos 6 ) 


u ,J •— u~ 


M 


2 Yp —(cos ff — cos 9 ) 

1 M + m ' ’ 


Also by the relation 


N .71 , T , 

1 + D + 1 + cl ~ N + n 

nd 


D = - 


N 4 (N 4 n) d 


and making these substitutions 


Yf 


Also 


F/= Nn [~J (^7 (l + D)* (1 + df 


7 Mm n „ N + n 

M 4 w " * P * N 4 (N 4 A) rf 

N 4 » 


(cos ff — cos 0) d 


, Ff . Mm , , r 
°g Yf ~ h M+ffl P N4(I4»)<I 


(cos 9 — cos ff) d. 


In forming (Yf — Ff) log (Ff/Yf) we see that the last factor is squared, and so 
the product contains the factor d 9 . We may, therefore, now write 1 for (1 + I)) ;} and 
(1 + df, and also write (N + w)/N for (N + «)/{N + (N + n) d}, otherwise we 
should have terms in d?. 

Therefore 


(Ff - Ff) log 


Nw (“) l (“) l e- /i f ai + w)V2 + ^ p ^ 


Id # 


Mm 
M 4 m 


4YV 


(N 4 nf 
~ N 2 


(cos 9 ' — cos 9 )' 2 . 


Again, dK/dt contains the factor \jj denoting the frequency of collision. Now the 
number of collisions given V and p in unit of volume and time is proportional to 
vdp, where s is the sum of the radii of M and m. Also all directions of p before 
collision are equally probable, and given the direction before collision, all directions 
after collision are equally probable. Therefore, given Y and p, the number in unit of 
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volume and time of colliding pairs for which the angle between Y and p before 
collision is between 0 and 0 d0 is Itts 9j p sin 9 d9, and the number for which after 
collision it lies between 6' and 9’ -j- d$" is jyirtdp sin 9' d9\ Hence we have to 
multiply (cos & — cos &f by rs 2 p sin 0 d0 sin 0' d&, and integrate in each case from 
7 r to 0. The result is § ns^p. 

And so we get 


dK 

dt 




Mm \ a /AMY* fhm\i 

M + m) \ IT ) \ 7T 


fco r it r2ir f 00 

X V 2 dV sin a da d/3 dp 4irp 9 

JoJoJo J o 


C Mm 


4Vy. 


77YT 


where V, a, /3 are usual spherical coordinates 


Also, as we have seen 
and therefore 


I dK 
K dt 


= — 0 suppose. 

And if K 0 , D 0 , d 0 be initial values, 


4~ £ (N + n 2 ) — !:) r, % d\ 

7T N ' ' (M + <sjh 

K = ||(N + n)d», 


16 (N +n) v/(M "' ) *** 


(M + m) 5 li 


K = K 0 e “ «, D = D 0 e d = d 0 e ~ ict , 

the rate of subsidence is directly proportional to the density and to the square root of 
the absolute temperature.* 


# Since writing the above I find that this result has already been obtained for the case of elastic 
spheres by Professor Tatt, by an independent method. 








